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The k=0 Friedmann Lemaitre Robertson Walker model with a positive cosmolog- 
■ ical constant and a massless scalar field is analyzed in detail. If one uses the scalar 

^ , field as relational time, new features arise already in the Hamiltonian framework of 

QJ \ classical general relativity: In a finite interval of relational time, the universe ex- 

Q ' pands out to infinite proper time and zero matter density. In the deparameterized 

quantum theory, the true Hamiltonian now fails to be essentially self-adjoint both 
in the Wheeler DeWitt (WDW ) approach and in LQC. Irrespective of the choice 
O of the self-adjoint extension, the big bang singularity persists in the WDW theory 

I I while it is resolved and replaced by a big bounce in loop quantum cosmology (LQC). 

j^jQl Furthermore, the quantum evolution is surprisingly insensitive to the choice of the 

self-adjoint extension. This may be a special case of an yet to be discovered gen- 
eral property of a certain class of symmetric operators that fail to be essentially 
, self-adjoint. 

p 

^ ; I. INTRODUCTION 

Loop quantum cosmology (LQC) of the k=0, A = Friedmann Lemaitre Robertson 
Walker (FLRW) model with a massless scalar field was discussed in detail in [1] . The scalar 
^ ! field serves as a viable internal time variable both in the classical and the quantum theory, 
^ I with respect to which relational observables such as the matter density and curvature evolve 
[2, 3]. This makes it possible to explicitly construct the physical Hilbert space and introduce 
relational Dirac observables to unravel physics of the Planck regime in a large number of 
cosmological models [4], and a scheme has been sketched even for full general relativity [5]. 
Using this setup it was rigorously established that, while the big bang singularity persists 
in the WDW theory of the k=0, A = model, it is resolved due to the quantum geometry 
effects of loop quantum gravity (LQG) [1].^ 
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^ Recently, this result has been conceptually sharpened using the consistent histories framework in which one 
can calculate probabilities for the occurrence of certain histories without recourse to external measurement 
devices or interaction with environment. Using appropriate coarse grained histories which completely 
decohere, it was shown that the probability of encountering a singularity in the distant past or future is 
1 in the WDW theory and in LQC for any state (which is in the domain of operators used to construct 
coarse-grained histories) [4, 6, 7]. 
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An appendix in [1] also outlined how the cosmological constant A with either sign can 
be incorporated. A subsequent, detailed discussion of the A < case appeared in [8]. It 
firmly established that, as in the A = case, in LQC the big bang singularity is replaced 
by a quantum bounce which occurs when the total energy density ptot reaches its maximum 
value pmax- Furthermore the numerical value of pmax is the same as in the A = case, 
Pmax ~ 0.41ppi, although now ptot includes a contribution from the cosmological constant in 
addition to the matter density p; ptot = p+A/SttG. It turns out that, by a suitable choice of 
time variable (or lapse function) , the A = model can be solved exactly [9] . This is not the 
case for A 7^ 0. Therefore, results of [8] for the A < case are conceptually important also 
because they demonstrate that the LQC bounce and the qualitative features of the resulting 
Planck scale physics are not tied to exact solvability. Finally, although the situation with 
the bounce is the same, the presence of the cosmological constant does alter the underlying 
mathematical structure in non-trivial ways. In particular, in the deparameterized picture, 
while the spectrum of the true Hamiltonian is continuous in the A = case, it is purely 
discrete in the A < case. 

The goal of this paper is to present an analogous, detailed account of the A > case. 
Even though we will again consider a massless scalar field, rather surprisingly, the flip of 
the sign of the cosmological constant changes the underlying mathematical and conceptual 
structure significantly. Let us begin with the classical theory. If one again uses the scalar 
field for internal time, in contrast to the A = and A < cases [1, 8], the Hamiltonian 
vector field on the phase space is now incomplete. As a result, volume of any compact 
co-moving region becomes infinite and the matter density vanishes at a finite instant 0o 
of internal time (j). This situation is qualitatively similar to that in the case of a non- 
relativistic particle in a steep negative potential whose dynamical trajectories reach infinity 
in a finite time. In such situations, typically, the Hamiltonian operator in Schrodinger 
quantum mechanics is symmetric but not essentially self-adjoint. Each self-adjoint extension 
then yields a unitary evolution but evolutions obtained from distinct operators are both 
mathematically and physically inequivalent. In the present case, one again finds that the 
true Hamiltonian operator generating evolution in the scalar field time is symmetric but 
not essentially self-adjoint. However, rather surprisingly, this ambiguity has negligible effect 
on states of physical interest: those that start out being peaked at a classical solution 
in a low curvature region. In particular, all these states undergo a quantum bounce and 
the total density ptot at the bounce is again universal. Furthermore, while the evolution of 
expectation values of physical observables does depend on the choice of self- adjoint extension, 
the dependence is extremely weak. This robustness may be related to the fact that, on the 
classical phase space, one can extend both the evolution equations and the solutions simply 
by analytical continuation, without having to introduce specific boundary conditions at 
infinity. Our analysis raises the possibility that there may well be a general pattern and 
new results could be found on properties of certain sub-classes of operators that fail to be 
essentially self-adjoint. 

The paper is organized as follows. In Sec. II we discuss the Hamiltonian framework for 
the k=0, A > model. Sec. HI is devoted to the WDW quantum theory and Sec. IV 
to LQC. We conclude in V with a brief summary and discussion. Because the numerical 
simulations in this paper were completed soon after the initial analysis in [1], they use an 
older value of LQC area gap which turned out to be half the value that is relevant for states 
used in LQC [10]. In the main text we use this more recent value. Therefore, unfortunately, 
in Sec. IV there is an occasional mismatch of factors of two between the text and the figures. 
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II. HAMILTONIAN FRAMEWORK 

Because our primary focus is on LQC, in this section we will summarize the phase space 
formulation of the FLRW model under consideration in terms of variables that descend 
from LQG. (For details, see, e.g. [4]). The space-time manifolds A4 will be taken to be 
topologically M^, equipped with a preferred foliation by spatially homogeneous and isotropic 
space-like slices M. The space-time metrics g will have the form 

g ^ -N'^dt^ + a'^(t)q'' (2.1) 

where q° is the fixed, positive definite, flat metric on M (determined by the co-moving 
coordinates), is the lapse function and a{t) the scale factor. Since all physical fields in 

the model arc spatially homogeneous, and since M is non-compact, integrals representing 
the symplcctic structure and Hamiltonians trivially diverge. One therefore introduces an 
infrared regulator — a cell C taken to be cubical with sides along co-moving coordinates — 
and restricts all integrals to it. The phase space structure and intermediate results depend 
on the choice of C and q°. Therefore one has to either show that the fi^nal physical results 
are independent of these choices or remove the infrared regulator by letting C to expand out 
to fill M. 

In LQG, one uses triads in place of 3-metrics. The freedom in the choice of their orien- 
tation enables one to introduce a configuration variable v which captures both the volume 
V of the cell C determined by q and the orientation of the underlying physical triad: 

where 7 is the Barbero-Immirzi parameter of LQG, = 47r\/37^pi the 'area gap' of LQG 
that is relevant to LQC [10] and Vg the volume of C with respect to q°. (sgnv is positive 
if the physical triad is oriented along the fiducial triad compatible with q° and negative if 
the orientations are opposite.) The canonically conjugate momentum is denoted by b. On 
classical solutions, it is given by 

b = ^XH = 7A-^ (2.3) 
' ' adt ^ ^ 

where H is the Hubble parameter and t is the proper (or cosmological) time. In these 
definitions, factors involving 7, A and £pi arc introduced to simplify the final expression of 
the Hamiltonian constraint operator later on.^ For the scalar field, the basic canonical pair 
is, as usual, 0,P(</,). Thus the phase space is topologically and equipped with basic Poisson 
brackets: 

2 

{b, v} = -, and {0,P(<^)} = 1 . (2.4) 

Because of a gauge fixing tailored to spatial homogeneity and isotropy, the Gauss and the 
diffeomorphism constraints of LQG are automatically satisfied [11]. Thus, we are left only 



^ The variable v is the same as that used in [1] and is related to the variable v of [4, 9] via: 1/ = \v and the 
variable b is related to b used in [4, 9] via b = (1/A)6. These relative rescalings by A make the variables 
u, b used here dimensionless which renders a considerable algebraic simplification in various expressions 
and equations. 
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with the Hamiltonian constraint. Since we want to use the scalar field (j) as internal time 
and since satisfies the wave equation on {M.^g)^ it is appropriate to use the lapse field 
that is adapted to a harmonic time coordinate r, satisfying Dr = [4]. As in [9], this is 
achieved by setting N :— a^. Then the Hamiltonian constraint is given by: 

C = - Snh^Gbh^ + TT'j^X^ h^G Av^ ^ 0, (2.5) 

Note that while b e (—00,00), on the constraint surface it must satisfy \b\ > bg — ^X^^/A/S. 
The equations of motion for the scalar field, generated by the Hamiltonian constraint, are: 

^ = and ^^ = 2p„„ (2.6) 

Eq. (2.6) implies that p(<^) is a constant in any solution and 0(r) = 2p(0)r + const. Therefore, 
in any space-time defined by a phase space dynamical trajectory, can be used as an 
evolution parameter in place of r. This is our relational time variable. The equation of 
motion for v in the relational time is given by 



[c^.P = 127rG[.^ + "''^P^ .1. (2.7) 



and its solution, expressed in terms of the physical volume V = {27TjXipi) \v\ of the cell C 
is: 

where 0o is a constant (that can vary from solution to solution). Eq. (2.8) implies that 
for each 0o we have two types of solutions: Those that start at = —00 with a big bang 
singularity, i.e. with zero volume for the cell C, and expand out to infinite volume at 
(f) — (po, and those that start out at infinite volume at = 0o and contract into a big crunch 
singularity at = 00 where C shrinks to zero volume. In either case, the evolution ends 
at the finite value 0o of the relational time 0. Put differently, the Hamiltonian vector field 
generating evolution in (f) is incomplete. 

An obvious question then is: Can the phase space evolution be naturally extended beyond 
(p — 00? We will now show that a mathematically natural extension does exist. Consider 
first matter density p = p'^^^/2V'^ of the scalar field, which is a physical observable of direct 
physical interest. Its time dependence is given by 

p(0) = ^ sinh2[x/T2^(0 - 0,)] . (2.9) 

Since it is analytic in 0, in the p-0 plane the dynamical trajectory represented by the 

'contracting branch' is simply an analytical continuation of 'expanding branch'. In space- 
time terms the full phase space trajectory can be interpreted as follows: the universe starts 
out with a big bang at (p = —00, expands out till the matter density p becomes zero at 
(p — (po, and then starts contracting, ending in a big crunch at = 00. Situation is similar 
with respect to b. The dynamical trajectory in the 6-0 plane also extends analytically from 
the expanding to the contracting branch: 



6(0) = ± 60 cosh[Vl27rG'(0 - 0o)] where bo := 7A y^ijs . (2.10) 
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Indeed, in both these cases, it seems artificial to stop the evohition of the expanding branch 
at = 0o and say that the contracting branch, (f) > 0„, is a distinct trajectory. 

In the space-time picture, on the other hand, as (p approaches 0o, the proper time goes 
to oo in the expanding branch and — oo in the contracting branch. Therefore the space- 
time represented by the expanding branch is future complete and that represented by the 
contracting branch is past complete. In effect, the extended space-time can be obtained by 
gluing together future null infinity of the expanding branch with the past null infinity of 
the contacting branch. Prom the perspective of space-time geometry, this gluing is just an 
optional mathematical construct. Prom the perspective of the Hamiltonian framework based 
on the relational time, on the other hand, the extension is natural and even appears to be 
necessary to have a complete picture of evolution. Since the Hamiltonian framework can be 
regarded as the imprint left on the classical regime by the quantum theory, one may suspect 
that the extension may have its true origin in the mathematical and conceptual framework 
underlying quantum cosmology. In the next two sections we will see that this is indeed the 
case. 



In this section we will first introduce the kinematical structure of the WDW quantum 
theory, then investigate properties of the operator ©a representing the gravitational part 
of the Hamiltonian constraint and finally discuss dynamics. As mentioned in section I, 
because of the presence of a positive cosmological constant, ©a fails to be essentially self- 
adjoint. Much of the discussion is devoted to establishing this property and exploring its 
consequences. 



In the Dirac program of quantization of constrained systems, one first ignores the con- 
straints and constructs a kinematical Hilbert space "Hkinj the quantum analog of the full 
phase space of the classical theory. The quantum constraint is then written as an operator 
on Hkin- Physical states lie in the kernel of this operator. 

In the WDW theory, one takes the kinematic Hilbert space to be T^j^jj^ = ® H^, 
where, as is common in the LQC literature [4], under-bars emphasize that the symbols 
refer to the WDW theory. As in the textbook Schrodinger quantum mechanics, one sets 
— I/^(R, d0) and :— L|(R, dw). Here the subscript S denotes that the states are 
symmetric, i.e., satisfy = '^{—v), thereby encoding the fact that since v ^ —v results 
from an orientation fiip of the physical triad , it is a large gauge transformation [3] under 
which physics of the model does not change. For later purposes we note that the inner 
product between two states ip, x the gravitational Hilbert space % is given simply by 



III. THE WHEELER-DEWITT THEORY 



A. Quantum kinematics 




(3.1) 



It is straightforward to write down the quantum operator corresponding to the constraint 
function (2.5): 
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C = I(8)9j + e^(8)I with (3.2) 
©A :=©o-7rG7^AWl and 0^ ^ -12nG^/\v\d^\v\dy^/\v\ , (3.3) 



where in the expression of 9^ we have chosen a factor ordering that is compatible with the 
one used in LQC in section IV (see [1] for details). 

Physical states lie in the kernel of C. To find this kernel and endow it with a Hilbert 
space structure, it is natural to use the general 'group averaging method' [12-14] as in 
[3, 15]. The implementation requires that the operator 0^ be self-adjoint and uses its 
spectral decomposition. In the A < cases these steps could be readily carried out and the 
physical sector of the theory could be constructed in a rather straightforward manner [1, 8]. 
With a positive cosmological constant, on the other hand, it turns out that the eigenfunctions 
ip{v) of Qj^ (with real eigenvalues) are Bessel functions of imaginary order. Unfortunately, 
for our purposes, they are difficult to work with (sec for example [16], for their properties). 
Therefore, it turns out to be more convenient to pass to the dual, 6-representation and work 
with the wave- functions ip{b). 

Let T> denote the dense domain in L^(M, dv) consisting of smooth functions which, to- 
gether with all their derivatives, fall off faster than any polynomial at infinity (i.e. let V be 
the Schwartz space). The operator is symmetric on V. Given any state ip{v) e V, the 
corresponding wave function in the b representation can be obtained by a Fourier transform. 
We will set 

[E^]{b) = I dv \v\-^ i^{v) e^"^ (3.4) 

so that, using the fact that vil^ip) = —2idbip{b) in the b representation, the operator 0^ 
assumes a convenient form: 

e^^-12nG[{bd,r-bld!]. (3.5) 

where bg — 7A ^/AJ3. We will first discuss certain properties of 0^^ and then use them to 
construct the physical sector of the theory. 



B. Properties of 9a 

1. Weak solutions to the eigenvalue equation 

To discuss dynamics, we have to extend 0a to a self-adjoint operator. To analyze exis- 
tence and uniqueness of these extensions, we need weak solutions to the eigenvalue problem 
— CV^) distributions such that 

(V^J©i-aix) = o VxeP, (3.6) 

where, as usual, the action (■0^|x) of the distribution (■0^[ on the test function \x) is defined 
using the Hilbert space inner product on "Hgr. 

Now, because of the factor |f |^/^ — introduced to simplify the form of the constraint 
operator — the form of the inner product in the b representation is not transparent. To spell 
it out, let us first note that the action of 0a leaves subspaces containing wave functions x 
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with support on positive and negative I'-axis separately invariant. Therefore, we divide the 
solutions to (3.6) into parts ip"^ 

t^^i^ = e{±v)±^ (3.7) 

by considering the inner product of V'c with test functions which have support only on the 
positive or negative parts of the f-line. Since by (3.1) the states supported on M+ are always 
orthogonal to the ones supported on M~ the action of can be written as 

(V^Jx) = (V^^lx)+ + (^^lx)_, where 

:= 4 / d6g(6) [±id,]x{h) Vx e I) . ^^'^^ 

With this explicit form of the action of at hand, we can now solve (3.6). For this, it 
is convenient to introduce a change of variables to simplify the form of 0^. Let us set 

^ rarctan(6/v^F^), |6| < 6„ 

\ sgn(6)[7r/2 - ln(6<,) + ln(|6| + ^/W^% 1^1 > ^o- 

or, reciprocally 

{ho sinx, \x\ < 7r/2, 

w, , . ^ / (3-10) 

OoCosh(|x| - -) sgn(a;), |a;| > 7r/2. 

where, as before, bo '■= 7A a/A/3. (Recall from section 11 that, in the classical solutions, 
\v\ ^ 00 as h ^ ±bo, or equivalently, as a; ^ ±7r/2.) Then, except at points b = ±bo we 
have^ 

= -127rG' sgn(|x| - 7r/2) . (3.11) 

Since Qj^ is just proportional to the simple operator except at x — ±7r/2, and we are 
interested in distributional eigenf unctions which are symmetric in x, they are necessarily of 
the form 

{A+e*^^ + A_e''^\ x>l 
B{e^^ + 6-^''), |x| < f (3.12) 



for constants A±, B± that satisfy suitable 'gluing conditions'. 

To determine these conditions we use the fact that the components ■^^ are independent 

solutions to (3.6) and apply the decomposition (3.12) of them directly to (3.6), using the form 
(3.8) of the inner product. Splitting the domain of integration of (3.8) onto three intervals 
Xj e { ] — 00, — 7r/2], [— 7r/2, 7r/2], [7r/2,oo[ } and integrating the resulting expression twice 



^ While 6a preserves the space of smooth functions of 6, on functions of x its action is discontinuous at 
X = ±7r/2 because, although 6 is a smooth function of x, db/dx = there. Since we are looking for 
distributional solutions, this discontinuity is harmless. 
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by parts we obtain 

se±i,(Te{+ -} ^ ' ' 



T^SinG Y] liin sgn{\x\ - n/2)[d^iljf]{x)[d^x]{x) 

se±i,<Te{+,-} 



±4:tJ2 f i9.x) [©A - CIV:^ , (3.13) 

where x {sit/2) denotes the hmit as x approaches sn/2 from above if o" = + and 
below if (T = — . Now, the integrand of the third term on the right hand side vanishes 
identically because -0^ are given by (3.12), and the second term on the right side also 

vanishes because smoothness (in b) of x implies that d^x at x = ±7r/2. Therefore, only 
nontrivial contributions to the right side of (3.13) come from the first term. Since Q/^x 
does not generically vanish at a; = ±7r/2 we conclude that weak solutions to the eigenvalue 
problem are given by (3.12) where the coefficients are chosen so that = iJj~^ + iJj~ is 

continuous in x = ±7r/2 (but not necessarily differentiable) . 



2. Self-adjoint extensions of Qj^ 

Deficiency Spaces: The operator is symmetric on T-L^ and the operator is 
known to be essentially self-adjoint [17]. However, the cosmological constant term acts like 
a negative unbounded potential. Therefore, from one's experience with Hamiltonians in non- 
relativistic quantum mechanics, one would not expect the operator 6^^ to be essentially self- 
adjoint. Indeed, its LQC counterpart was recently shown to admit a family of inequivalent 
extensions [18]. We will now show that this general expectation is correct by analyzing 
deficiency spaces [19] of 0^. For notational simplicity, we will first rescale this operator and 
consider 9^ = (127rG')-i 9^. 

The deficiency spaces JC^ are spanned by (kinematically) normalizable solutions -0^ to 
equation (3.6) with the eigenvalue ( — ±Si ^. Their elements are solutions to (3.11) which 
are symmetric and everywhere continuous. Therefore, from (3.12) it follows that they take 
the general form 

f^{x) = Ae2(iT^)l"l + 5e-2(i^^)l^l, (3.14) 

for \x\ > tt/2, where A,B are some constants. However, we will now show that only the 
solutions with A = are normalizable. Heuristically this is plausible because the first term 
on the right hand side diverges as |x| — >■ oo while the second goes to zero in this limit. But 
to establish the result we need to consider the Hilbert space norms. Since the inner product 
is simple in the v representation, let us solve the equation directly in that representation. 



The deficiency spaces as defined in [19] correspond to tire eigenvalues ( = zLi, however one can equally 
work with the spaces corresponding to ^ = ±ir where r is any positive real number. We chose r = 8 just 
to simplify notation in subsequent calculations 
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Using from (3.2) we can express ijj as a linear combination of Bessel functions 

^^{v)^C ^ '1- ^ ^' +D ^ 'L_ ' (3.15) 



Since these functions decay at infinity sufficiently fast, to find if they are normalizable it 
suffices to focus on their behavior for small v. In this limit, the solution approaches 

-2(l=Fi)ln|«| 2(l=Fi)lnb| 
7P^{v)^E_ ^^+E+ (3.16) 



whence from the form of the inner product (3.1) it follows that only the solution with E_ = 
is normalizable. Since this normalizable solution goes as |f |^^^, its Fourier transform (3.4) is 
well-defined and square integrable in b. Having established this property, we can now return 
to the eigenfunctions (3.14) in the x representation. Since db/dx ~ sinh(x), the normalizable 
solution iIj{x) must have the property that J da;sinh(a;) < oo. This condition implies 
that the normalizable solution is given by (3.14) with A = 0. 

Setting A — and using the symmetry properties of elements of 'Hgi and the continuity 
property of eigenfunctions we conclude that elements of deficiency subspaces /C^ have the 
form 

{l\p2{l±i)x I „-2(l±i)a;l |„| ^ tt 
' . X f.n M , ^ " ' (3-17) 

cosh(7r) e-(2(iT0kl— )^ k| > f ^ ^ 

for some B & C. From this it follows immediately that /C^ are 1-dimensional. 

Self-Adjoint Extensions and their Domains: The fact that K,^ are non-empty 
immediately implies that admits a family of inequivalent self-adjoint extensions [19]. 
Elements of this family are labeled by the unitary transformations U : /C"*' — > IC~ . In our 
case they are all of the form 

Ua-.ip^ ^ e^>-, (3.18) 

where are some chosen normalized elements of /C^. Thus, 0'^ admits a 1-parameter 
family of self-adjoint extensions, labeled by a G [0,27r). Using theorem X.2 of [19], we 
conclude that the domains Va of these extensions are given by 

= {^ + ^+ + P, 3/^^ e /C^}. (3.19) 

Up to a constant rescaling, the terms ijja '■— ijj~^+U°'ijj~^ € "Da that depend on and characterize 
the extension are given by 

. ..^Ul/2)[e'^cos{2x-a/2) + e-'^cos{2x + a/2)] , \x\<n/2, 
-"^ ^ |cosh(7r)e-(2N-'^)cos(2|x| -q;/2), \x\>7r/2. ^' ' 

Although these functions are continuous, generically, they are not differentiable x — ±7r/2. 
On the other hand, each element '0 of I? is differentiable in b and satisfies [dxil^\{x — ±7r/2) = 
0. Therefore, for fixed value of a the ratio between the left and right hand derivative of any 
element of Pq, is a constant, common for all depending only on a. Furthermore, the 
function 
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(where denote derivatives from above and below respectively) is monotonic: 

= °t9?'°"r\ ^ (3.22) 
cosli(27r) — sin(Q;) 

Therefore the relation 

[0,27r) 9 Q;^/3(a) e [0,7r) (3.23) 

is a bijection. Thus ^ provides an alternate labeling of the extensions which is much more 
convenient as it has a direct interpretation in terms of the discontinuity in the derivatives 
^^■0 at a; = ±7r/2. Therefore, labeling the extension by ^ allows us to easily identify the 
eigenfunctions of which span particular extended domain. 

To summarize, the symmetric operator 0^ defined on the domain T> admits a 1-parameter 
family of self-adjoint extensions ©a,/?) parameterized by ^ e [0,7r) with domains given by 
(3.21). The parameter [3 directly captures the boundary conditions at a; = pi/2 satisfied by 
states in the domain of O^^^ via (3.21). 

Eigenfunctions: Given a self-adjoint extension ©a,^ we can construct the corresponding 
physical Hilbert space V^^^ by group averaging and discuss quantum dynamics. Both these 
tasks require us to find the eigenfunctions of ©a,^- For this, we first note that, since the 
operator d"^ in (3.2) is negative definite, only the positive part of ©a,^ relevant for the 
construction of ■ The general, symmetric eigenspaces of ©a corresponding to the eigen- 
value are 2-dimensional and are spanned by functions (where A; e M) of the general 
form 

(^) = A UosHkx), \x\ < 7r/2, 

^ \cosh()t7r/2)e^^(l-l-'^/2), \x\ > 7r/2 ^ ^ ^ 

where A & C The eigenfunctions spanning the particular extended domain need to 
satisfy the condition analogous to (3.21). They are therefore given by 

-^'''^ ^ \cos{k\x\+ap{k)), \x\>7r/2, ^ ' 

where B e ]R+, k e IR"^ and the phase shift (Jp{k) satisfies the relation 

iMk^/2 + a,{k))^'-^^^. (3.26) 

To find the normalization factor B wc use an analog of the method used in appendix A2 
of [20]. Specifically, first note that, for large \x\, ip^^ approaches an eigenfunction of 0a=o 
with a controlled rate of convergence 

±l^(x) - B cos[fc(ln(2) -ln(M + | + 5) + ^^(A;) ]. ^^'^^^ 
where x = In Similarly, in the ^-representation, we have [16] 

^/3,.(^)=^;,.(^) + 0(^')- (3-28) 
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which imphes [20] that and satisfy the same (distributional) normahzation condi- 

— p,re — p,K 

tions. Therefore, the known normahzation of "0° [9] and the relation between -0° (x) and 

— p,k — p,k 

. ("y) determined by (3.4), (3.9) fixes B as 

B = -^=. (3.29) 

To summarize, for every self-adjoint extension ©a,^ of 0^ the positive part of its spectrum 
Sp+(0A/3) equals the entire positive half of the real line and is absolutely continuous 
and non-degenerate. The domain of the operator is spanned by eigenfunctions ^, 
with k e M"*", given by (3.25), where the coefficient B is fixed by (3.29). The resulting 
eigenfunctions satisfy the normalization condition 

iep,k'k0,k)=Hk-k'). (3.30) 



C. The physical sector 

Knowing the self-adjoint extensions of can trivially construct all the self-adjoint 

extensions of C_ by substituting (3-2) with ©a,/? = 127rG'0A_^. Then, knowing the 

spectral decomposition of ©a ^ one can find the physical Hilbert space corresponding to 
each extension by the group averaging [12-14] method as in [3]. The resulting (positive 
frequency) physical states corresponding to each extension are of the form 

*(x,0)=/ dfci(fc)e^,,(x)e-«^ (3.31) 

Jr+ 

where cu(k) = \/V2tiG k. As explained in [3, 9] the physical inner product given by the group 
averaging procedure reduces to simply 

(*il*2)=/ d'^4i(^)i2(^) so that T^P^y = L2(R+,dA;). (3.32) 

Knowing this form of physical states, we can introduce as usual [4, 9] the notion of the 
evolution with respect to the internal time 0, and interpret the Hamiltonian constraint as 
the evolution equation. The resulting unitary evolution of the initial data at = 0o is given 

by 

^(x, (t)o) ^ *(x, 0) = e'^^-'^°^ *(x, (Po) , (3.33) 

where |0a,/3| is the positive part of ©a,^- 

To extract physics from this setup, one needs a family of relational observablcs on liphy, 
parameterized by the internal time 0. In the cases A < [1, 8] this role was served by 
operators l^l,^ corresponding to the volume (of the fiducial cell C) at given value of 0. In 
our case, however, in the classical theory v becomes infinite at a finite value of (f). On the 
quantum level this property is reflected by the fact that an operator \v\^ can map even the 
elements of the Schwartz space S G H^^^ outside of "Hp'^^. Since the physical origin of this 
technical problem is clear from the classical analysis, so is the solution: One can simply 
replace \v\^ with |/('fi)|</,, where / is a bounded function of v. A convenient choice of / is an 
'angular parameter' 9 e]0,7r/2[ defined by 



(3.34) 
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where K is a constant with dimension of the vohime. Following the general procedure 
introduced in [1, 3], one first introduces a kinematical operator acting on the initial data 
space P^iLgr consisting of states in 2igr in the positive part of the spectrum of ©a,/3 

{9±) (v) = (P+ arctan(i^) P+ ^) {v) (3.35) 

{P^ projects elements ip{v) of "Hgr into their positive frequency part by restricting the 
corresponding Fourier transform t/j{k) to the positive half /c-line.) In the second step, one 
extends it to the physical operator ^|<^ using (3.33): 

(^U„ M/) (,T, 0) = e^^^-"^") (e^) (x, 0,) , (3.36) 

This relational Dirac observable enables us to effectively track the evolution of volume of 
the fiducial cell C. 

There are other observables which are manifestly independent of the choice of a fiducial 
cell. An example is provided by the family \b\^ built out of the kinematical observable 
This observable is of direct physical interest because classically b is proportional to the 
Hubble rate H: b = 7A H. Finally, another cell-independent observable of direct physical 
interest is the total energy density ptotU commonly used in LQC [9, 15]. It is constructed 
starting from the kinematical operator 

We will now use this setup to explore quantum dynamics of the WDW theory. 



D. Quantum dynamics 

Thanks to the simple form of the evolution operator (3.11) and its eigenfunctions (3.25), 
we can draw some general conclusions regarding quantum dynamics at a semi-heuristic level. 
Note first that, for |a;| > n/2. every self-adjoint extension 0a,/3 yields a standard Klein- 
Gordon equation in the a;, plane, with boundary conditions at the 'barrier' \x\ = 7r/2 
defined by that particular choice of (3. Most of the wave is refiected at this barrier as the 
tunneling amplitude is exponentially suppressed. Therefore one would expect the qualitative 
features of the evolution to be as in Fig.l. Consider a quantum state ^(k) which is sharply 
peaked about a large k* for which "^{x, (p) is peaked on a classical trajectory 

6(0) = bo cosh[Vl27rG'(0 - 0^)] or a;(0) = ^ + VUnG {(t)o - 0) (3.38) 

in the low curvature region at some < 0^. Then, because of the refiection at = ±7r/2, 
one would expect ^(a;, 0) to remain sharply peaked on the trajectory 

x(0) = ^ + Vl27rG|0-0„|. (3.39) 

for all 0, even to the future of = 0o. But this is only a qualitative argument because 
it neglects the tunneling into the region \x\ < 7i/2 and ignores the k dependence of the 
reflection coefficients. 



13 




X 



FIG. 1: Qualitative behavior of the quantum evolution can be inferred from the classical dynamics 
coupled with the form of eigenfunctions. Expanding classical trajectories start at (x = c«, </> = —oo) 
(the big bang) and end at {x = tt / 2, cj) = (j)o) (where the matter density p goes to zero), while 
the contracting trajectories start at (x = it/2, (p = 4>o) and at (x = cxd, </> = cxd) (the big crunch). 
For definiteness we have set (/>o = in the hgure. The form of 0a,/3 ^'^d its eigenfunctions in 
the x-representation suggest that incoming semi-classical wave functions would be largely reflected 
at X = 7r/2 and follow the contracting classical trajectory which is the analytical continuation of 
the original one. The figure shows only the part x > because the physical states ^!.(x, 0) are 
symmetric under reflection in x. 



To make precise statements regarding the global evolution one has to analyze the expec- 
tation values and dispersions of appropriate observables. Since the model is not analytically 
soluble, one has to resort to numerical methods using a suitable class of states. Following 
the existing literature, we select the Gaussians sharply peaked about some large k* and some 
h* at given initial time 0o <^ —\j\fG. The spectral profiles corresponding to such states are 
of the type 

|r(fc) = e"^^^ ^ ^here 0* := 0o + ^ arcosh(676o) . (3.40) 

V 127rG' 

Since the explicit form of the basis functions e^j^. is known via (3.25) and (3.29) the wave 
function ^(6, 0) can be calculated by first carrying out a direct integration of (3.31) and 
then passing to the h representation using (3.9)). 

Recall that the expression of the inner product is quite complicated in the 6-representation 
but extremely simple in the ^-representation. On the other hand, the expression of some of 
the basic operators is simpler in the 6-representation. Therefore the calculation requires us 
to pass between the two representations using the Fourier transform (3.4) and its inverse. 
Consequently, numerical analysis of dynamics — calculation of the expectation values and 
dispersions of ^|</,,ptotU — was carried out in the following steps 

(i) First, the profiles ^(x, 0) at a given were found via the integration of (3.31) over 
the interval G [fc* — 7(T, fc* + Tcr] via Romberg method with the stepsize and order 
selected dynamically to achieve the relative integration precision of the order 10~^. 
The corresponding wave function ^(6,0) was then calculated using (3.10). 
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FIG. 2: Quantum dynamics in the WDW theory. Evolution of the wave packet initially peaked 
at = 5 X 10^ and |^|^=o ~ """/^ with dispersion Ap^/p^ = 0.03. Figures (a) and (b) show the 
evolution in the v-cj) and h-cj) plane, respectively. Figures (c) and (d) show the expectation values 
of the observables \9\^., and /StotU in this state and compare them with the analytically extended 
classical trajectories. In these simulations, /3 = 0, the constant K in the definition of 9 was set to 
is: = 5 X 10^, and A = 



(ii) Next, ^(6, (j)) was used to compute the function 

$(6,0) := 62-^(6,0). (3.41) 

that is needed to evaluate the expectation values in step (iv) below. 

(iii) Results were transformed to the f- representation via an inverse of (3.4). The Fourier 
transform involved in it was computed via a Fast Fourier Transform (FFT) algorithm 
with 2^ probing points (where N varied from 16 to 19). To adapt the resolution in 
h (corresponding to the size of the domain of calculation in v) the domain of trans- 
formation h G [— &Af , was allowed to change in the process of evolution and varied 
from 0.5 to 5. 
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(iv) Finally the norm, and the expectation values were calculated as follows 

11*11^ = / dv|*(v,0)|^ (3.42a) 

(I^U) = ll^ir' / arctan(M)|^(^;,0)|2^ (3.42b) 
Jr 

(ptotU) = Cp\m~^ [ dv^iv,cf))^v,(l)), (3.42c) 
where Cp — 3/(87rG A^7^). The dispersions were found via a standard relation 

Si = id') - {of (3.43) 

for a given observable O. On the other hand, the expectation values of squared oper- 
ators are determined via 

(1^15) = ll^ir' / arctan2(M)|^(^,0)|2^ (3.44a) 
Jr 

{pl) = C'p\m-' f dvmv,<P)\'. (3.44b) 

All the integrations in this step were performed via the trapezoid method, with the 
configuration of the probing points preset by the FFT in the previous step. 

A large number of simulations were performed varying the parameters of the initial state 
and the self- adjoint extension used in evolution. They bore out the semi- heuristic expecta- 
tion described in the beginning of this sub-section. An example of the results is presented 
in Fig. 2. For any choice of extension, states that are sharply peaked in the distant past (or, 
future) remain so throughout the evolution and follow the analytically extended classical 
trajectories towards big bang and big crunch singularities. Thus, as in the A = case, the 
big bang and the big crunch singularities are not resolved in the WDW theory. The new 
element is that the wave packets follow the analytically extended classical trajectory beyond 
(f> — (f)o, irrespective of the choice of the self-adjoint extension.^ 

Now, usually in quantum mechanics the existence of incquivalcnt self- adjoint extensions 
implies that at some point the physical evolution breaks down and additional data are 
required to continue it. In our case these data correspond to the gluing conditions at 
\b\ = bo or |a;| = n/2. However, surprisingly, for states which start out as semi-classical in 
the weak curvature regime, this non-uniqueness does not manifest itself at least at the level 
of the expectation values of the observables: All of them follow the (analytically extended) 
classical trajectories dictated by parameters k*,b*, (j)* used in their construction. This result 
is a reflection at the quantum level of the fact that the classical solutions admit a unique 
analytic extension. 

But what about generic wave functions? To better understand the dependence of physical 
quantities on the choice of extension in this case, let us first consider wave functions in the 



^ A previous result in the WDW theory, but using a 'spinor formahsm' [21], hinted at the existence of such 
an extension. Although the scale factor played the role of time there, the quantum wave packet developed 
a second branch which, in retrospect, may be interpreted as the analog of this extension. 
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X variable. For |,x| > tt /2 each basis function ^, has the form of a standing wave and can 
be decomposed onto the incoming e^^ and outgoing e^^ component 

e% = {l/VMk) e±^(^l^l+<^Mfe)). (3.45) 

By replacing e^^ in (3.31) with e^^, we arrive at a split of the wave function '^{x^ 0) onto 
expanding and contracting components. Let us define 

eH^) = ^""'^'^'^ e%,{x) (3.46) 

by rotating each basis function so that the result §!.^{x) is /9-independent and rewrite the 
expanding and contracting wave functions in terms of them: 

•^(x, (/))=[ dA;i'^(fc) e;^{x) e^'^(*^)('^-'^°) . (3.47) 
Jm.+ 

Then the spectral profiles satisfy 

|r'^(A;) = e^'''^^'''>±{k) , whence i'^(A;) = e^'"^^''^^' {k) . (3.48) 

Thus upon reflection, the initial expanding wave proflle ^~{k) undergoes a phase shift 
g2!o-^(fe) ^i^Q^ depends on the choice of the self-adjoint extension. Note that the expression 
(3.26) of (Tp implies that the difference — ct/j/ in rotations of the phase caused by two 
different extensions is not global, but depends on k. This subtlety should be reflected in the 
evolution of observables. 

To extract this information, it is convenient to follow [20] and regard the global evolution 
as a scattering process. For, in the distant past, i.e., (f) <S —1/VG, and in the distant future, 
i.e. (f) ^ 1/\/G, we are near the big bang and big crunch singularities where the effect 
of the cosmological constant can be neglected relative to that of matter, i.e., A/SnG ^ p. 
Therefore, in these 'asymptotic' regions, dynamics is better and better approximated by 
that in the A = case. Thus we can regard 'incoming' and 'outgoing states' as belonging 
to the A = physical Hilbert spaces. The incoming state can be thought of as being 
scattered because of the presence of the cosmological constant which dominates dynamics 
near |a;| = tt/2. As in [20], the scattering process can be analyzed using observables In \v\^. 
However, it is more convenient to consider observables In deflned in analogy to \b\^. The 
limit of the basis functions is provided already by (3.27). Following our analysis in the x 
variable, we split this limit into incoming and outgoing components expressed in terms of 
the basis of in the 6-representation 

k{b) = ^ ei(fc/3o+a,(fe)) g+(^) ^ ^-i{k0o+a,ik)) g-^(^) ^ (3(^-2^ ^ ^3 49^-) 

v2 

with efe(6) = (l/V47r|A;|)e^'=i"l''l, (3.49b) 

where /3o := 7i/2 + ln(2/6o). Using this split in (3.31) we can deflne the asymptotic future 
and past states ^ of the spectral profiles: 

^^{k) = e{^k) e±^(l*^l/3o+<^MI*^l)) ±{\k\) . (3.50) 
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Now the global evolution is described by the scattering operator S 

1^') ^ 1:*"^) = -51^-) (3.51) 

where the matrix elements of the scattering operator S are given by 

S{k, k') := (e;.| S le^O = g-^^^g-C^OC^o+'^^dfcl)) ^(^ + A;'). (3.52) 

The extension dependence is encoded entirely in the phases shift crp{k) given by (3.26). For 
large A;, it can be expanded as 

<^m--^ + l-fi-\ sin(2/3)e-^- + 0{e-^'-) . (3.53) 

Let us now select two extensions corresponding to some /3 and /3' and two states \^), \^) 
corresponding to those extensions, sharply peaked at some k* and such that in the asymptotic 
past they are equal. Then by (3.51) in the asymptotic future they are related just by a phase 
rotation 

|r+(A;) = ^-2i{af,{k)-a^,{k)) |r'+(;t). (3.54) 

This allows us to estimate via (3.53) the leading order correction to the difference between 
the expectation values of In as 

1(^+1 ln|6|^|:*+)-(^'+| ln|6|^|:*'+)| =7rcos(/3 + /3')sin(/3-/3')e-'^'=*+ 0(e-2"'=*). (3.55) 

(Since and are asymptotic states, i.e. states in the A = theory, the difference 
is 0- independent.) To appreciate the physical implication of this result, let us focus for 
a moment on semi-classical states peaked at a large k* such that initially the dispersions 
in ln|6| and ln|A;| are comparable and the uncertainty product is approximately saturated. 
Then, in either the /3 or /3' theory, because of the Heisenberg uncertainty principle, the 
dispersion in In \b\ is of the order of \k*\~^^^^\ Because of the exponential suppression in 
(3.55), the difference in the expectation values of In \b\ in the /3 and theories is therefore 
completely negligible compared to the dispersion in In |6| in either of these theories. 

Let us now return to general states. The dependence of dispersions of the choice of 
self-adjoint extension is slightly stronger than that in (3.55). Results of [20] show that 
the possible growth of dispersion between the asymptotic past and future states depends 
on the behavior of dkCTp{k). By repeating for the WDW theory the derivation of certain 
triangle inequalities on dispersions obtained in [20] for LQC, one can shows that the expected 
difference between the dispersions 5+ and 5^ of In |6| for the two states \^)-, \^^) (of the 
WDW theory now under consideration) will be of the order of 

^\[dl{ap-ap,)]{k^)\a,^\k\ (3.56) 
< ^ I cos(/3 + (5') sin(/3 - (5') \ e''^^^* 5,^\^. 

where 5\n\k\ is the dispersion of the observable In |A;|. (Here and in what follows the approx- 
imate sign in the inequality emphasizes the fact that we are keeping track only the leading 
order terms.) 
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Now, the change in the dispersion between the asymptotic past and future state is given 

by 

\5+-5^\<2{^+\A(dkaf,) |*+) 

<2|[9|a^](r)Uin|fe| (3.57) 
<7r2sin(2^) e--'=*5in|fe|, 

where S- is the dispersion of the observables ln|6|^ on the state |^~). (Again since the 
asymptotic states refer to the A = theory, this difference is 4> independent.) Thus, the 
difference in the dispersions — in the asymptotic future in the /3 and (3' theory is 
comparable to the change in the dispersion — between asymptotic past and future 
in any one theory. Both these quantities are very small compared to the dispersion 5in|jt| in 
the scalar field momentum. 

To summarize, computer simulations of states which start out as Gaussian in the weak 
curvature region established that the big bang and the big crunch singularities fail to be 
resolved in the WDW theory. They also brought out the surprising fact that although the 
final physical sector of the quantum theory docs depend on the choice of the self- adjoint 
extension 6a,;3 of the symmetric operator 0^, the difference in the dynamics of these states 
is neghgible. The S-matrix strategy first introduced in [20] enabled us to obtain certain 
analytical inequalities for generic states for which the dispersion in k — i.e., in the field 
momentum — is finite. They showed that if the dispersion 5in|fc| is small, then differences in 
the S-matrix predictions of physical theories that result from different choices of self-adjoint 
extensions are enormously suppressed. 

IV. LOOP QUANTUM COSMOLOGY 

In this section we will show that the procedure followed in Sec. Ill for the WDW theory 
can be repeated in a rather straightforward manner for LQC. Therefore our discussion 
will be parallel to that of Sec. Ill A and HID with emphasis on the differences from the 
WDW theory. 

A. LQC kinematics 

As in the WDW theory, the kinematical Hilbert space is a tensor product Ti^^^ — "Hgr <S> 
"H^. However, while we again have = L^(R, d0), the gravitational Hilbert space is now 
different [4, 11]. As in the A = case [1], it is given by "Hgr = I/^(]R, d/^Bohr), where M is a 
Bohr compactification of the real line and d//Bohr the Haar measure thereon. A convenient 
basis is again provided by the eigenvectors of the operator i): 

v\v) = v\v), so that V\v) = (27r7A4i) l^"! 1^") , (4.1) 

where, as before, 7 is the Barbero-Immirzi parameter of LQG and is the LQC area 
gap [10]. As in the WDW theory, in the volume representation states in Tigr become wave 
functions il>{v), which are again taken to be symmetric iplv) = ip{—v) to incorporate the 
fact that V — > — V is a large gauge transformation corresponding to the flip of the orientation 
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of the physical triad. However, unhke in the WDW theory, the ip{v) now have support only 
on a countable set of points along the v-axis and their inner product is given by a sum 

(^1^0 = J2 i^Wiv) , (4.2) 



rather than an integral. 

Because the Hilbert space is so different, the differential operator — 
— 127rG -\/|?7[9y \v\ dj;y^\v\ of the WDW theory fails to be well-defined on T-Lgr- Therefore, 
one has to first express the classical constraint in terms of the elementary variables of LQG 
— holonomies and fluxes — and then promote the result to an operator on "Hgr. This sys- 
tematic procedure leads to the following form of the constraint operator [1] 

(7 = I (g) a| + Ga <8) I, Ga := eo-7^G'7VA^;^ (4.3) 

where 

-[eo^l^]{v) = U{v) i^{v - 4) - fo{v) i^{v) + f-{v)i^{v + 4) , (4.4) 
with the coefficients /o,± given by 

f±{v) = (3716-/4) ^v{v±A) {v±2), fo{v) = {2,1:0/2^ . (4.5) 

Thus, the second order differential operator of the WDW theory is now replaced by a 
second order difference operator Oo with uniform steps of size v = ±4. Therefore, there is 
super-selection: one can investigate dynamics separately on uniform lattices in the f-space 
and each sector consisting of wave functions with support on any one of these lattices is 
preserved by the complete set of Dirac observables of interest, discussed in section III. In 
this paper, we will restrict ourselves to the lattice C = {v = 4n, n e Z} for simplicity 
because in LQC physical results are largely insensitive to the choice of the sector [22]. 

Finally, as in the WDW theory, for technical reasons it is more convenient to work in 
the dual representation in which states are wave functions ip{b) of the conjugate variable b. 
However, since the LQC states ip{v) have support only at v — An, we now have a Fourier 
series in place of the Fourier integral (3.4): 

[^V^](6) = -^ J2 k|-^(^)en (4.6) 

where the point v = was removed from the transform because the state with support just 
at f = is dynamically decoupled from the orthogonal sub-space spanned by states which 
vanish at v — 0. Since are supported on Cq, their images J-'ip are periodic in b with the 
period tt. Therefore one can restrict the support of the wave functions [J'^](6) just to the 
circle b e [0, tt], with the identification [7^-0] (0) = [^V](^)- 

By inspection, the elementary operators v and jV"^ defined by 

v\v)=v\v), and J\f^\v) = \v + fj.) , (4.7) 
in the v representation are transformed to 

v^2idb, and Af^^e-''^/''. (4.8) 
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in the b representation. As a consequence, the operator ©a assumes the form 

Oa = -127rG' [ {sm{b)dtf - by,] , (4.9) 

in the b representation, where, as before bo := "yX a/A/3. 

Let us first consider the case when 60 > 1 or A > Ac := 3/7^A^. In this case, 0a 
is essentially self-adjoint, whence one can readily repeat the procedure of section IIIC to 
construct the physical Hilbert space. However, because this Oa is negative, the physical 
Hilbcrt space is now zero dimensional! (For proofs, see [18].) Thus, in striking contrast to 
the WDW theory, in LQC a non-trivial quantum theory exists only when the cosmological 
constant A is less than a critical value, A^. Although this result is not phenomenologically 
relevant because Ac is of Planck scale, it is of considerable conceptual interest. In the rest 
of this section, then, we will with work A < Ac. 



B. Properties of Ga 

1. Weak solutions to the eigenvalue equation 

Note that, in the 6-representation, the gravitational part of the Hamiltonian constraint 
is a differential operator as in the WDW theory. This suggests that one may be able to 
simplify it by a change of variables. We will now show that Ga can in fact be transformed 
to the same form as in the WDW theory: 

Oa = 1271^ sgn(|x| - Xo) dl . (4.10) 

However, there are two key differences. First, now the new variable takes values on a compact 
interval x e [—xm-iXm], with points —xm and xm identified. Second, the transformation is 
much more complicated in that x is defined in terms of the elliptic integral of the first kind 
F{y,k): 

' -j^^ Fib', 1/(1 -bl)), \b'\<Bo, 

\xM-lF{7r/2-b',l/bl), b'>Bo, (4-11) 
-XM +°f F{7i/2 + b', l/c2), b' < -Bo, 



where 



b' :— b — 7r/2, Bg :— arcsin(6o) 



=F(7r/2 - Bo, 1/(1 - bi)), XM ■■= Xo + -rF{Bo, 1/6^) . (4.12) 



The dependence of Xo and xm on bo is shown in Fig. 3. In particular 

lim Xo{bo) = +00, lim [xm — Xo]{bo) = n/2, (4.13a) 

]im Xo{bo) — n/2, lim [x m — Xo]{bo) — +00. (4.13b) 

bo-^l bo^l 

For the technical simplicity, it is convenient to introduce the rescaled variable 

y.^TT— e[-7r,7r], (4.14) 
Xm 
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(a) 



(b) 
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FIG. 3: (a) The dependence of the functions Xq and xi := xm — Xq (of (4.12)) on the cosmological 
constant A. (6) The dependence of Uo = t^Xq/xm on A. Here Ac = 8/7^ is the critical value of A 
at which the energy density in the cosmological constant equals the maximum energy density pmax 
allowed by LQC in this model. 



in terms of which the operator ©a takes the form 



G/ 



X 



Ga, where 6^ := sgn(|?/| - yo) d\ 



M 



(4.15) 



with Uo := ttXo/xm- 

We are now ready to analyze weak eigenfunction of 0'^. As in the WDW analysis 
we will have to go back and forth between the v and the y representations. The weak 
eigenfunctions are distributional solutions to 

(V-cie't - ci|x) = o Wxev , (4.16) 

where V now consists of states x ^ ^gr which have support only on a finite number of points 
of the lattice C on the v-axis. The particular form of the coefficients /o,± of Go in Eq. (4.5) 
and the form of the inner product of 'Hgr in the f-representation allows one to split each 
such eigenfunction -ip,^ into components ip^ 



^lJfiv) ■.= ei±v)^^iv), (4.17) 

which again satisfy (4.16) as in the WDW theory. Secondly, as a function of y, any solution 
to (4.16) is necessarily of the form 

A{e'^y + e 



A.), 



\y\ e [0,2/0) 
\y\ e (2/0, vr] 



(4.18) 



analogous to (3.12). These two observations allow us to directly apply the techniques devel- 
oped in the final part of Sec. Ill A to find the necessary and sufficient conditions for (4.18) 
to satisfy (4.16). The result is a complete analog of the one in the WDW theory: the weak 
eigenfunctions have to be continuous in y = ±yo but not necessarily differentiable there. 
This property transfers directly to the operator Ba for which any (weak) eigenfunction is 
again globally continuous, but not necessarily differentiable with respect to x at a; = ±Xo. 
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2. The self-adjoint extensions of 0a 

We know already from the analysis of [18] that 0a admits a one parameter family of self- 
adjoint extensions, labeled by elements of U{1). Each extension corresponds to particular 
asymptotic behavior of the basis functions e{y) in the limit v — >■ oo. Recall, however, 
that in the WDW theory, 6-representation allowed us to find a more useful interpretation 
of the choice of extension in terms of gluing conditions at 6 = ±6o. The similarity of the 
WDW and LQC constraint operators and form of the eigenfunctions when expressed in 
terms of x variables (introduced respectively via (3.9) and (4.11)) suggests that a similar 
interpretation should exist also in LQC. We will now show that this expectation is correct. 

For simplicity we again consider the operator 0'^. Its deficiency functions ip"^ G /C^ 
are the weak solutions to the equation (4.16) with the eigenfunctions ^ = ±i. Because of 
symmetry and the global continuity are of the form 

^^M = cl^'/^^)^^''^'''^ + ^"''^'''^)' f4 19) 

^ |(l/c_)(e(i±*)(-2')/^ + e-(i±*)(-^)/^) \y\>yo, ^ ' 

where 

c+ = e(i=F*)s'°/^ + e-(iT*)f°/^, (4.20a) 

C_ = e(l±i)(T-2/o)/V2 _^ g-(l±^)(7r-J/<,)/^/2_ (4.20b) 

Since all the eigenspaces are non-degenerate, there exists a 1-1 correspondence between 
these solutions and the appropriate eigenfunctions of 0a discussed in [18]. Since it was 
shown in [18] that all eigenfunctions of 0a are normalizable, it follows that eigenfunctions 
(4.19) are also normalizable. The unitary transformations between /C+ and /C^ and the 
extended domains are again given by the exact analogs of (3.18) and (3.19) respectively. 
The extension-characteristic terms V'a -analogs of (3.20)- are now of the form 

^.(x) = C'|^(^'^-")' (4.21) 
/(7r-|/,7r-7/o, -a), \y\> yo, 



where C" e C and 



fiy, yo, a) := ^ ^ V e(2/+<^2/.)p/v^ cos 

[cosh(x/2|/„) + cos(x/2|/„)] ^^^^ 



ayo _^ pa 



V2 2 



(4.22) 



Since [dxip]{x — ±Xo) = ^ip e V, the extensions are uniquely determined by the 
parameter (3, 

R( \ + ([dyMiy = yo)\ ^ . . . 

'''"^ ^= "^'"^ IWM^J 

that encode the ratios between the left and right derivative of ■^a at the gluing point y = yo- 
With this parametrization, the choice of self-adjoint extensions can be directly interpreted 
in terms of gluing condition at 6 = ±6o also in LQC. (Recall that in the classical theory, 
this is where the universe reaches the infinite volume.) Furthermore for any fixed value of 
yo e]0,7r[ the function (3{a) is a bijection of the circle of radius 1 onto the circle of radius 
1/2. (See Appendix A for a proof of this assertion). Therefore, as in the WDW case, /3 is 
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a convenient label for the extensions, equivalent to a, and Eq. (4.23) succinctly captures 
the domain of 6a,/3- These properties will be used in the next subsection to identify 
the spectra of particular extensions 6a,^ and the corresponding eigenbases in the physical 
Hilbert spaces. 



C. The LQC physical sectors. 
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FIG. 4: (a) The lowest three eigenvalue ujn are plotted as functions of the self-adjoint extension 
parameter (3 and compared with their approximations co'^ obtained by neglecting the remnant 
Q(-g-27m(7r-?/„)/?/o) (4.27), for A = 0.01 Ac. (b) the frequency gap Aoj is plotted as a function of 
cosmological constant. 

We can now fix any one self-adjoint extension Oa,/3 of Ga in the quantum Hamiltonian 
constraint (4.3) and perform group averaging to obtain the physical Hilbert space Ti^^ . 
Because the spectra of 6a,/3 are discrete [18], the physical states are of the form 

oo 

^(x, <p) = Y, e/3,„(x) e^-^-(<^-^°), (4.24) 

n=0 

where the spectral profiles are square-summable sequences and e^^n are the normalized 
eigenf unctions of 9a,/3 (belonging to P^), with eigenvalues Alternately, e^,„ are eigen- 

functions of the operator 6^^^ with eigenvalues k'^, related to ujn via Un = {VI^ttG tt/xm) kn- 
They satisfy Eq. (4.23) and their functional form is given by 

g (a;) = ]V Uosh[kn{n -Vo)] cos{kny), \y\ < Vo, ^^25) 
|cos[/c„?/o] cosh[A;.„(7r - y)], \y\ > i/o, 

where A^/3,n is a normalization factor and the eigenvalues A;^^„ are determined by the condition 

tan(A;„?/o) + tanh[A;„(7r — i/o)] tan(/3) = 0. (4.26) 
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The form of (4.26) implies, in particular, that the spectra of Oa,^ for different values of f3 
are disjoint. Furthermore the lowest non-negative eigenvalue q is an invertiblc function 

of f3, whence the physical Hilbert spaces corresponding to each extension are in fact 
different subspaces of "Hgr. 

Although the eigenvalues A;^ are provided only imphcitly, the relation (4.26) allows us to 
determine their asymptotic behavior for large n explicitly: 

kn = (nTT - f3)/yo + 0(e-2-"('^-2/o)/2/o)_ (4 27) 

Consequently, as n grows, the distribution of kn quickly approaches the uniform one with 
the separation Ak = tt/Uo- This property transfers directly to 0a, where 

hm [ojri+i - oJn] = Au; := VUnG n^/xo, (4.28) 

Although the distribution of lower frequencies is not quite uniform, it approaches uniformity 
extremely fast (see Fig. 4a), and the asymptotic separation Aw depends only on the value 
of the cosmological constant (see Fig. 4b). In particular 

hm Aw = and lim Aw = 27rVl27rG. (4.29) 

The form of the quantum Hamiltonian constraint (4.3) and of the physical states (4.24) 
allows us again to use the scalar field as an internal time and regard the square root of the 
positive part of ©a,/3 as the generator of time evolution on the subspace Ti^^ — P^Tigr of 

^{x, 4)o) ^ ^{x, 0) = e*('^-^°)\^®^ ^{x, 4)o) V *(x, 0) e P/?^gr, (4.30) 

where, as before, is a projection onto the positive part of the spectrum of ©a,/3- 

To obtain the physical consequences of this evolution, we will use observables |^|^ and 
PtotU in exact analogy with Eqs. (3.35) and (3.37) of the WDW theory. However, in contrast 
to the WDW theory, we cannot define \h\^ since the operator dy fails to exist in LQC; we 
can only define its bounded, periodic functions such as sin( 



D. LQC dynamics 

To facihtate comparison with the WDW theory, we will use states closely resembhng the 
Gaussians used in Sec. Ill D peaked about large w*: 

where 0* is given by (3.40). we can then repeat the procedure used in the WDW theory 
to carry out numerical simulations. Specifically, these computations were performed in the 
following steps: 

(i) First, the spectrum Sp(6a,/3) was found from (4.26). In all cases considered, the 
approximation provided by the analytical expression (4.27) turned out to be excellent. 
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(ii) Given any eigenvalue in Sp(Ga_/3) the corresponding cigcnfunction ipp^niv) was found 
in the domain v & {An : n — 1 . . . Ni} hy solving the iterative difference equation 

ujliPp^niv)^[eA,pip0,n]iv) (4.32) 

with initial data ip/3,n{v = 4) = (—1)". In various simulations, Ni ranged from 10^ 
to 2.5 X 10^. Using the expression (4.2) of the scalar product in the I'-representation, 
the norm of ■0/3,n('t') was calculated using a polynomial extrapolation. The normalized 
eigenfunctions e^,„(f ) were then used together with the profile coefficients ^„ of Eq. 
(4.31) to obtain the physical state ^(^,0) in the ^-representation: 

oo 

^{v, <P) = J2^n e^A'") e'^^'-^'^-^"^ (4.33) 

n=0 

where the sum extends over those n for which uj^ € [<^* — 7cr, u)* + la] . This state is 
defined on points v e {4n : n = 1 . . . A^i} and is simply the transform of the state in 
Eq. (4.24) to the I'-representation. A key difference between the WDW theory is that, 
once the eigenvalues are obtained using (4.26), we can work directly in the v repre- 
sentation because, unlike in the WDW theory, the LQC eigenfunctions e/3^„(w) in the 
V representation can be easily handled numerically. (Recall that in the WDW theory, 
the corresponding eigenfunctions are Bessel functions of imaginary order which are 
difficult to deal with numerically.) 

(iii) Expectation values and dispersions of the relational Dirac observablcs, |^|^, ptotU were 
computed directly in the v representation using the expressions (4.7) of v and the 
shift operator M^. In particular, when the operator ptot is factor ordered as in the 
WDW theory, it becomes a simple combination of the shift operators: 

Note that the expression has the same form as 3.37) in the WDW theory, but the 
WDW operator is now replaced by the bounded operator sin(6) . 

As in the WDW theory, a large number of simulations were performed by varying param- 
eters of the initial state and the self-adjoint extension used in quantum dynamics. Fig. 5 
illustrates the results of a typical simulation. The qualitative behavior is the same as that 
in the A < case [8]. 

• States under consideration remain sharply peaked over a large number of 'epochs', 
where each 'epoch' is characterized as the evolution between consecutive quantum 
bounces. Thus, in each epoch the universe starts out with the maximum but finite 
total density pmax, expands out till the matter density vanishes and then recoUapses, 
the density again reaching pmax at the end of the epoch. 

• Expectation values of the Dirac observables \d\^ and ptot\(j> are well approximated by 
the classical effective dynamics discussed in Appendix B throughout the evolution, 
including the quantum bounces. The difference between the two is much smaller than 
the dispersions of the corresponding Dirac observables. 
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FIG. 5: LQC dynamics. Figures (a) and (b) show the evolution of the (absolute value of the) 
physical wave function in the v-(p and b-cp planes respectively. Figures (c) and (d) compare the 
expectation values of the Dirac observables \6\^ and ptotU with those in the effective and classical 
theories. Away from the Planck regime there is excellent agreement with (the analytical extensions 
of) solutions in general relativity. But in the Planck regime (left ends of figures (c) and (d)) there is 
a very large departure from the classical behavior because of quantum geometry effects. However, 
the effective trajectories capture the quantum evolution very well even in the Planck regime. In 
these simulations, the initial state was a gaussian peaked at = 5 x 10^ and l^l^^g ~ arctan(lO) 
with the relative spread Ap^/p^ = 0.03. The self-adjoint extension corresponded to /3 = 0, the 
constant K in the definition of 9 was set to K = 5 x 10^ and the cosmological constant was 
A = O.OlAc ^ 0.2£pf . 



• For A <^ Ac := 3/7^A^, quantum dynamics is well-approximated by the (analytically 
extended) classical general relativity trajectory so long as the scalar field energy density 
is small compared to the Planck scale, p|<^ ^ ppi. 

• However, as p|tot = pU + A/(87rG') approaches the Planck scale, an effective repulsive 
force due to quantum geometry starts to dominate. It soon overwhelms the classical 
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gravitational attraction and forces a quantum bounce at ptot — Pmax- As a result, the 
contracting quantum universe bounces back into an expanding quantum universe. As 
noted above, the wave functions remain sharply peaked even during the bounce and 
expectation values of Dirac observables follow those given by the effective equations 
even in the Planck regime. 

• As in the WDW theory, quantum evolution is not unique because it depends on the 
choice of the self-adjoint extension ©a,/3 of ©a- Each extensions corresponding to the 
choice of a specific refiective boundary condition, but at ptot — ^/i^'^G), and not in 
the Planck regime near ptot ~ Pmax- 

• For any chosen self-adjoint extension, as in the WDW theory, the expanding universe 
recollapses after reaching zero matter density. However, unlike in the WDW theory, 
the big bang and the big crunch is replaced by quantum bounces. Therefore, we now 
have a nearly periodic evolution in LQC admitting an infinite chain of the bounces at 
the Planck energy density and recollapses at zero matter density The duration of each 
'epoch' is given by 

to an excellent approximation. 

• However, the distribution of cj„ is not exactly uniform. Consequently, quantum states 
slowly disperse from one cycle to the next. But for large a;„ deviations from uniformity 
decay much faster those in the case of A < considered in [8]. Therefore the rate of 
dispersing is much slower than the already low rate found in [8] . 

• As in the WDW theory, the dependence on the choice of the self-adjoint extension is 
surprisingly weak for the states considered here. For generic states, results of the S- 
matrix discussion of [18] show that the dependence on /3 of expectation values of Dirac 
observables is negligible compared to the dispersion of the corresponding quantities. 
The dependence of dispersions on (3 is also very small compared to the dispersion 6\n\k\ 
in the scalar field momentum (which is a constant of motion). 



V. DISCUSSION 

In this paper, we analyzed in detail the WDW theory and LQC of the k=0, A > FLRW 
model along the lines of the treatment of the A < case of [8], thereby completing the 
program outlined in an Appendix of [1] . As in the A < cases, the scalar field can be used 
as a global clock, providing us with a natural notion of relational time both in the classical 
and quantum theories. However, interestingly, there is a key difference in the physically 
most interesting case, that of A > 0: In classical general relativity, solutions that start with 
infinite matter density p at the big bang at time 4> — —oo expand out and now achieve p = 
at some finite value 0o of internal time (when the volume v of the fiducial cell C becomes 
infinite). Thus, in the p-0 plane each of these dynamical trajectories starts out at (f) — — oo 
but ends at = ^o- But it can be analytically extended beyond (f) = (po and the extension 
represents a universe which starts out with zero matter density at = 0o but contracts, 
ending in a big crunch singularity at = cxo. ^From the relational time perspective, then, 
one is led to regard the two branches as providing a single dynamical trajectory because 
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it is artificial to simply end dynamics at a finite value of time. Now, in non-relativistic 
mechanics if the potential is negative and steep the particle may roll off to infinity in a finite 
amount of time. In that case, one has to choose from a one (or more) parameter family of 
boundary conditions at infinity to continue dynamics beyond that time. In the present case, 
by contrast, we did not have to resort to making a choice because the most interesting Dirac 
observable, p|^, is analytic in (p. 

Nonetheless, the fact that the universe expands out to p = at a finite value of the 
relational time introduces ambiguities in the quantum evolution: The operator ©a which 
generates dynamics with respect to now fails to be essentially self-adjoint. Before dis- 
cussing this point in detail, let us first note two aspects of this phenomenon. First, it is 
not a peculiarity of LQC; it occurs also in the WDW theory. Second, in both cases, the 
lack of essential self-adjointness is related directly with the behavior of the system at large 
V and low matter density p; its origin does not lie in the Planck scale physics. In both 
quantum theories, the dynamical operator admits a one parameter family of self-adjoint ex- 
tensions. For a general system, different choices of extensions can give rise to very different 
dynamics. However, in this model the results are surprisingly robust with respect to this 
choice. Not only is the qualitative behavior of dynamics the same, but the differences in 
the dynamics of the expectation values of the most interesting Dirac observables in theories 
resulting from two different extensions are smaller than their dispersions in any one theory, 
even for general states. Furthermore, numerical simulations show that, irrespective of the 
choice of extension, quantum states which are semi-classical in the low curvature (or low 
total density ptot) regime remain sharply peaked at the eodended classical trajectory in the 
low curvature regime both in the WDW theory and LQC. It is tempting to conjecture that 
this robustness of quantum dynamics is related to the fact that we did not have to choose 
a boundary condition at = 0o to extend the classical p — 4> trajectory. In the remainder 
of this section, most of our discussion on the behavior of wave functions will refer to these 
states. 

As in the A < cases, there is a pronounced difference between the quantum dynamics 
of the two theories in the Planck regime. In the WDW theory, the wave function simply 
follows the extended classical trajectory into the big-bang and the big-crunch singularities. In 
LQC, by contrast, while these states remain peaked at the classical trajectory so long as the 
curvature (or ptot) is low compared to the Planck scale, there is a dramatic departure in the 
Planck regime. There is again a new repulsive force with origin in the quantum geometry that 
overwhelms classical gravity and cases a quantum bounce. Again, the numerical simulations 
show that, although the force is so strong in the Planck regime, it dies very quickly and 
becomes negligible once ptot falls below 10~^ — 10~^ Planck density. In LQC then, even 
though we are in the k=0 case, we are led to a scenario that is approximately cyclic. As in 
the k=l LQC models, the quantum evolution spans an infinite number of epochs. In each 
epoch the universe begins with a quantum bounce where ptot ~ 0.41ppi, expands out till 
Ptot = A/SttG and then undergoes a collapse till it reaches another quantum bounce. For 
states under consideration, dynamics is nearly periodic. 

How does this dynamics appear in the space-time picture? Let us begin with the classical 
theory and consider a solution in which the universe starts out with a big-bang at = — oo. 
It expands out to J+ — which is space-like for A > — where the matter density vanishes 
and Ptot = A/SttG. In terms of the physical metric, this space-time is future complete. 
However, the extended phase space trajectory analytically continues the space-time geometry 
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across effectively gluing it with J~ of a contracting solution.^ Quantum states under 
consideration remain peaked at these extended space-time geometries across J. An extension 
is but to be expected both in the WDW theory and in LQC: quantum evolution in the 
internal time is unitary and (p achieves a finite value (po at J of the given classical solution, 
unitary evolution could not just stop there. What is interesting is that, irrespective of the 
choice of the self-adjoint extension, the state remains sharply peaked on the analytically 
extended geometry. 

This extension, and the ensuing nearly cyclic scenario has some similarities with Penrose's 
recent proposal of a cyclic conformal cosmology [24]. However, there are also key differences. 
In our case, J"*" of the expanding branch is glued to the 3~ of the contracting branch; not 
to the big-bang singularity of the next 'aeon'. More importantly, quantum geometry effects 
are crucial in LQC. In particular h appears in the denominator of the expression of the 
maximum density pmax whence, as one would expect, pmax would diverge in the classical 
limit h ^ 0. Therefore, quantum effects and a non-zero h play an essential role in the 
approximately cyclic scenario of LQC. By contrast, a central feature of the cyclic conformal 
cosmology paradigm is that, although one does have unboundedly large curvatures, h plays 
no role at all in this regime. 

Use of the scalar field as a relational time variable played a key role throughout our 
analysis, both in the classical and quantum theory. What would have happened if we had 
made some other choice? In a recent analysis [25, 26] non-rotating dust has been used in 
place of the scalar field. In this case, the expression (4.3) of the gravitational part Ga of 
the Hamiltonian constraint is modified because the lapse is now tailored to proper time. In 
particular, the coefficient of A is now linear rather than quadratic in v. Consequently, the 
analog of Ga is now essentially self-adjoint and the LQC evolution resembles that in the 
A = case [1]: The universe starts out with infinite volume in the distant past, collapses, 
undergoes a quantum bounce and then expands out to infinite volume. However, in the 
Planck regime, quantum matter should be described using quantum field theory and for 
all standard quantum fields the kinetic term in the Hamiltonian is quadratic in momenta. 
Therefore the specific feature that simplifies the mathematics in the case of dust is no longer 
available and the overall situation is then the same as that in the case of the scalar field. 

But what if we return to using scalar field as matter source but let a geometric variable 
be the relational time? An obvious choice is volume. But in LQC, volume fails to be single 
valued making it difficult to introduce the 'time-dependent' relational Dirac observables 
— such as the matter density operator p\y — especially near the bounce. On the other 
hand, the conjugate variable b is single valued both on classical and effective trajectories 
and, as the form of the Hamiltonian constraint suggests, the function y it determines is a 
possible candidate for relational time. However, it appears that the evolution would then 
be unambiguously unitary only for \y\ > y^. Furthermore, defining the physical state would 
require specification of the initial data aX, y — yo but the theory does not provide any 
selection principle for this task. These difficulties with v and h could well be surmountable 
with new ideas and more careful analysis. But as of now they seem to be more serious 
handicaps than the complications associated with the use of the scalar field as relational 
time we encountered in this paper. 



The detailed gluing procedure will involve a conformal completion along the lines of [23] where the normal 
component to J of the metric is rescaled by a different power of the conformal factor than the tangential 
one. 
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Could one perhaps retain scalar field as the matter source but set lapse = 1 in the 
quantum constraint and use a 'timeless framework' in LQC? Results of [15, 27] imply that 
we would have been led to a theory that is mathematically free of ambiguities associated 
with self- adjoint extensions. However, as we now explain, this theory is difficult to interpret 
and it is unclear whether it is physically viable. We will conclude our discussion with a 
detailed elaboration of this point. 

With lapse = 1, the quantum Hamiltonian constraint has the form C = B{v) (8> i9| + 

® I where Cgr is the gravitational part of the Hamiltonian constraint and B{v) ~ 1/v 
for large values of v (see, e.g., [1]). This total constraint operator C has been shown to 
be essentially self-adjoint [18] and for simplicity we will denote its self-adjoint extension 
also by C. One can therefore use it directly for group averaging and construct the physical 
Hilbert space "W^^ in the 'timeless framework' without recourse to deparametrization. Then, 
although states can be represented as wave functions ^(v, 0), the physical scalar product is 
no longer given by an integral over t; at a fixed value of 0. What is the relation between this 
W^'^ and the Hilbert spaces 1-1^^^ associated with self-adjoint extensions ©a,/3 we constructed 
in this paper? It turns out [27] that "H^^^ is huge; it is given by the direct integral of all 
^Phy. ^phy _ ^^^-Hfy^li, where / is the interval (0,7r). 

We will briefly discuss a simple example — due to Wojciech Kamiiiski in [28] — to illustrate 
the relation between these Hilbert spaces. Consider a 2-dimensional strip, M = M x [0, 1], 
with coordinates e R and x e [0, 1] and a constraint thereon in the form of the Schrodinger 
equation {—id^ + idj:)'^{x, 0) = (with a first-order Hamiltonian). As in LQC, plays the 
role of time while x is to be thought of as the analog of the compactified volume coordinate 
9. The operator '■= —idx fails to be essentially self-adjoint on the closed interval [0, 1]; 
it admits a l-parameter family of self-adjoint extensions Qx,p-, labeled by /3 G [0, 2n) (with 
domain given by wave functions ^^{x) satisfying ■0(1) = e*^'0(O)). For each extension, we 
can construct a physical Hilbert space T-lF^^ in the standard manner from solutions to the 
quantum constraint: ^!{x,(f)) e H^^^ if and only if —id^'i>{x,(f)) — ©j;,^ ^(x, 0). Next, let 
us define "H := J^jTif^^ df3, with / = (0,27r). ("H is analogous to Ti^^^ obtained by group 
averaging in [27]). Every ^^(x, 0) e yJg"^ can be expanded as 



oo 




(5.1) 



n=— oo 



where k := 27rn -|- /3. On the other hand ^(x, 0) G "H has the form 




(5.2) 



where we have set k — 27rn -|- ^ as above. The norms in T-L^^^ are given by: 




(5.3) 



Note that "^{k) has support on the entire /c-axis and furthermore, the norm of "^{x, 0) also 
involves the integral of |^(/c)P over the entire k axis. Therefore, while elements ^/3(x,0) 
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of any one T-L^ ^ are restricted to have support only on the physical configuration space 
M = M X [0,1], in Ti^^^ are allowed to be non-zero all along the ,x-axis, even 

though points outside the x-interval [0, 1] have no physical interpretation in the model under 
consideration. In other words, although we restricted ourselves only to a 'strip' M = Mx [0, 1] 
of the Minkowski space in defining the classical system, in effect W^^ describes a system 
on the entire 2-dimensional Minkowski space M^. ^ To summarize, for the particle on the 
strip M, one would physically expect that the quantum theory should be formulated just 
on M since values of x outside [0, 1] have no physical meaning. This expectation is borne 
out if one works with any one self-adjoint extension Qx,i3 of but not if one works with 
the direct integral W^^ of all the resulting Hilbert spaces. 

In LQC, the situation is analogous. Working with a specific self-adjoint extension allows 
us to remain in the interval [0, oo] of the |f |-axis, introduce Dirac observables and track 
their evolution in the internal time (p and compare it with classical trajectories. On the 
other hand, working in the timeless framework in effect requires us to extend the I'-axis 
beyond v — oo and this extension is difficult to interpret physically.^ For the same reason, 
while one can introduce Dirac observables also in the timeless framework, we cannot ask 
for their 'evolution' and it is difficult to compare predictions of the quantum theory with 
those of the classical. In particular, while the choice N — 1 oi the lapse yields evolution 
in proper time in the classical theory, unfortunately this interpretation does not extend to 
the quantum theory in a simple way. Thus, while at first it seems mathematically natural 
to work with the lapse = 1 because the full constraint is then essentially self-adjoint, the 
resulting Hilbert space l-L^^^^ appears to be simply too large to be physically viable in the 
above context. 

Finally, our analysis brought out an unexpected robustness of the quantum evolution 
with respect to the choice of self-adjoint extensions Qa,i3- As we noted above, this may be 
related to the fact that, in the physical sector associated with any self-adjoint extension, 
under unitary evolution quantum states of interest follow the natural and unambiguous 
analytic extension of classical trajectories. Is this perhaps a special case of as yet unknown 
general result? Is the lack of sensitivity of the quantum evolution on the choice of analytic 
extensions have its origin in some special features of the classical evolution? 
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Note that, since the basis functions e'(^'^"+''^^ can be extended analytically to the entire x-axis, the 
quantum constraint can be uniquely extended to the full Minkowski space K^. 
^ Furthermore, unlike in the simple example discussed above, here one cannot extend the constraint operator 
uniquely even if one expresses the basis vectors as functions of the 'compactified volume' coordinate 9 of 
(3.34), as these functions are not analytic in ^ = 7r/2. 
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Appendix A: Bijectivity of I3{a) 

In section IV B we began by labeling the self-adjoint extensions of 0a by a parameter 
a and then switched to a more convenient parameter (3. In this appendix, we will show 
that the map a — )■ ^ is bijective, i.e. that the (3 parametrization used to denote self-adjoint 
extensions as ^a,/? is viable. 

The definition of (3 and the periodicity /3{a) = (3{a + 7r) — a direct consequence of its form 
given by (4.21), (4.22) and (4.23) — imply that f3{a) is a well defined continuous function 
mapping from a circle of radius 1 to a circle of radius 1/2. Therefore to establish the desired 
bijective property it suffices to prove that, for any given i/o e]0, 7r[, the derivative of (3 with 
respect to a is bounded and isolated from zero. 

Let us consider first the function X{yo, a) := tan(/3) which can be decomposed as follows 

X{yo,a) = r I I 7, (Al) 



where 



/i = cosh(V2(7r - y,)) + cos(V2(7r - 7/„)), (A2a) 
/2 = cosh ( \/2yo) + cos ( \/2yo) , ( A2b) 



and 



Qi = sinh(A/2yo)[cos(Q;/2) — sin(Q;/2)] 
— sm (y2yo)[cos(a/2) + sin(a/2)], 

g2 = sinh(y2(7r - yo))[cos{a/2) + sin(a/2)] 
+ sm{a/2- V2{n -yo)). 



(A3a) 
(A3b) 



It is straightforward to check by inspection, that both fi and /3 are strictly positive and 
isolated from 0. 

The derivative of X over a takes the form 

^ yyo, a = r X 2/ T' ^4 



where the function 



/s = cosh(y2(7r - 2yo)) - cosh(y27r) 

+ cos(^(7r - yo)) sm{V2yo) - smh{V2yo) (A5) 
+ sin(-\/2(7r - yo)) sm{\/2yo) + smh.{V2yo) 



has (also by inspection) the following properties 

• h{yo = 0) = 0, 

. Vy„e [10-3,7r] Myo) G [-50,-0.1], 

. Vy„ e [0,10-3] [dyJs]{yo) G [-123.0,-122.5], 
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which again implies its boundedness and isolation from zero for any < i/o < n. 
The derivative d^P can be expressed using X, X' as follows: 

and by the properties of fi, f2, fs it never vanishes and can reach infinity only when both 
gi and g2 vanish simultaneously. The condition gi = g2 — can be reexpressed as the 
linear system of equations for sin(Q;/2) and cos(a/2) (treated as independent variables). 
The existence of the nontrivial solution to that system requires that the determinant of the 
equation matrix vanishes, that is 

[sinh(v^(7r - yo)) + cos(v^(7r - y^))] x [sinh(V^yo) - sm{V2yo)] 

+[sinh(v^y„) + sin(v^y„)] x [sinh(v^(7r - - sin(^(7r - = 0. 

Since all the terms in the square brackets are explicitly positive within G]0,7r[ this con- 
dition is never satisfied. Therefore, the denominator of the rightmost expression in (A6) 
never vanishes, so for fixed value of yo the derivative da/3 is bounded due to the continuity 
of X, X' and the compactness of the domain of a. Thus we have 

> c > Miyo, a)>C> -00. (A8) 



Appendix B: Classical effective dynamics 

In Sec. IV D we found that wave functions which start out being sharply peaked on a 
classical trajectory in the low curvature region remain sharply peaked throughout the evo- 
lution, including the Planck regime. This strongly suggests that there may be effective 
dynamical trajectories on the classical phase space which incorporate the appropriate quan- 
tum corrections in the Planck regime and approximate the full quantum evolution quite well 
throughout evolution. A particularly convenient heuristic method to arrive at these effective 
equations was proposed in [29] and later derived more systematically analytically in the case 
A = [30]. It was successfully tested in several cases of isotropic LQC [1, 8, 9, 31]. Here we 
brieffy recall how these effective equations arise by adapting the detailed discussion for the 
A < model from [8] to the A > case now under consideration. 

The strategy comes from a geometrical formulation of quantum mechanics in which the 
quantum Hilbert space is regarded as an infinite dimensional symplectic manifold, sometimes 
called 'the quantum phase space'. The idea is to find an embedding of the finite dimensional 
classical phase space into this infinite dimensional 'quantum phase space' such that the full 
quantum evolution preserves the image of the embedding to a good approximation. For a 
harmonic oscillator, the embedding is provided by coherent states (whose dispersions are 
determined by the mass and the spring constant). Such an embedding is possible also for 
FLRW models [30, 32] (for a brief summary, see [4]). The result is a quantum corrected, 
effective constraint. Although it is simply a function on the classical phase space — obtained 
by taking the expectation values of the quantum constraint operator in states corresponding 
to the image of the embedding — it differs from the classical constraint function by terms 
involving h. As mentioned above, dynamical trajectories generated by this effective Hamil- 
tonian constraint have turned out to provide an excellent approximation to the full quantum 
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evolution of states that start out to be sharply peaked around a classical trajectory in the 
low curvature regime. 

In the LQC literature, the effective Hamiltonian constraint is written using lapse = 1 
so that it generates evolution in proper or cosmic time. Therefore, to facilitate comparison, 
we will do the same here. Then, the effective Hamiltonian constraint is given by 



n 



off 



\p\2 sin (jjc) + 



+ 



3 

p2 



8nG-f^jj,^ ~"' ' 2\p\l ' 8nG 



A^O. 



(Bl) 



By calculating the Hamilton's equations from Ties, and using (Bl) again to simplify the 
resulting expressions, one arrives to the following evolution equation for the energy density 
p of the scalar field: 

1/2 



p' = ±4^3^ 



PPtot 



1 - 



Ptot 



Pn 



(B2) 



where ptot '■= P + A/(87rG) and the prime denotes the derivative with respect to the scalar 
field 0. (Thus, one first calculates the derivatives of p and with respect to proper time by 
taking Poisson brackets of p and with "Hgff and then combines them to find p' and then 
p'.) However this equation is inconvenient to use in numerical simulations because it is not 
regular at p = and ptot = Pmax- Also, the sign in front of the righthand side changes in the 
process of evolution (at the bounce where ptot = Pmax and at the recollapse where p = 0). 
Therefore for the purpose of finding the solution it is more convenient to use the second 
order equation derived from (B2), 



p" = 2A'kG 



(2p + Ptot) 1 - 



Ptot 
pc 



p 



(B3) 



which admits a unique global solution to the initial value problem with the initial data 
{p{(t>o), p'{4>o)), where p'{(t)o) is determined from ptot via (B2). 



1.6 



(B6) 
(4.35) 



1.2 1 



0.4 



FIG. 6: The period A0 of the evolution predicted by the effective theory via (B4) is compared 
against the approximate period of the genuine quantum evolution (4.35). 

Although this equation can be solved analytically, its solution is expressed in terms of 
the elliptic integrals of the first and second kind whose values have to be found numerically. 
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Therefore, it is more convenient to solve it numerically from the beginning. In the actual 
calculations (such as the ones presented in Fig. 5) this step was carried out with the use of 
the adaptive Runge-Kutta method of the 5th order known as RK45 (Cash-Carp). 

The form of equation (B3) implies already that the trajectory p(0) is periodic, with the 
period given by a direct integration of (B2): 

^ rptot=Pc 

\/V2ttG Jp=o 

The dependence of A0 on the value of cosmological constant is presented in Fig. 6, where its 
values were calculated numerically via the standard trapezoid method. The extremal points 
correspond respectively to the values p = (minimum corresponding to the recoUapse at 
infinite volume) and p — Pmax — A/ (SttG) (maximum corresponding to the quantum bounce). 



PtotP 



1 - 



Ptot 



(B4) 
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